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The emergent universe scenario provides a possible alternative to bouncing cosmology to avoid 
the Big Bang singularity problem. In this paper we study the realization of the emergent universe 
scenario by making use of Quintom matter with an equation of state across the cosmological constant 
boundary. We will show explicitly the analytic and numerical solutions of emergent universe in two 
Quintom models, which are a phenomenological fluid and a nonconventional spinor field, respectively. 
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I. INTRODUCTION 

Inflation is considered as the most successful model 
of describing physics of very early universe, which has 
explained conceptual issues of the Big Bang cosmology 
[TJ-[3] (see [4-6 for early works). Among these remark- 
able achievements, inflation has predicted a nearly scale- 
invariant primordial power spectrum which was later 
verified in high precision by Cosmic Microwave Back- 
ground (CMB) observations [Jj. The success of infla- 
tion is mainly based on a series of assumptions includ- 
ing an enough long period of quasi-exponential expansion 
and the applicability of perturbation theory during this 
phase. However, these assumptions often bring troubles 
to inflation models, such as the fine-tuning problem of 
the potential parameters. Moreover, it was pointed out 
that inflation models suffer from the initial singularity 
problem inherited from the Big Bang cosmology [8]. 

Recently, there are increasing interests in alternatives 
to inflationary cosmology which can not only be the same 
successful as inflation in explaining early universe physics 
but also avoid the initial spacetime singularity (for exam- 
ple see [9] for a recent review) . One class of the alterna- 
tive scenario is the bouncing cosmology, which suggests 
the expansion of our universe was preceded by an ini- 
tial contraction and then a non- vanishing bouncing point 
happened to connect the contraction and the expansion 
[TUrfT^] . It was pointed out in Ref. [T3] that a real- 
istic cosmological model realizing a nonsingular bounce 
requires a matter field with Quintom behavior |14j of 
which the equation of state (EoS) of the universe has to 
cross over the cosmological constant boundary w = — 1 
twice. This type of bounce model was extensively studied 
in the literature in recent years (for example see [T5HT8] 
and references therein, and see |19j for a review on Quin- 
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torn cosmology). A bouncing model combining the Mat- 
ter Bounce scenario |201 121] and Quintom scenario |14j 
was achieved in the frame of Lee- Wick cosmology |22j . 
which showed that such kind of bounces can give rise to 
a scale-invariant power spectrum for primordial curva- 
ture perturbation. Later, the cosmological perturbation 
theory of bouncing cosmology was established to a com- 
plete frame which includes primordial non-gaussianities 
[2"51 [24] , entropy fluctuations and curvaton scenario [25J , 
and related preheating process |26j . A recent nonsingu- 
lar bouncing model |12j which combined the benefit of 
Matter Bounce [2Ql [21] and Ekpyrotic cosmology [27] [28] 
was proposed and the w crossing —1 scenario can be re- 
alized without pathologies through the Galileon-like La- 
grangian (see for the Galileon model and [301 EI] f° r 
extensions) . 

Besides the bouncing, there is another interesting cos- 
mological scenario which could be a strong competitor 
of inflation, which is the so-called emergent universe sce- 
nario [32-34 . This scenario suggests that our universe 
was initially emergent from a non- vanishing minimal ra- 
dius and experienced a enough long quasi-Minkowski 
phase and then entered the normal thermal expansion. 
It was obtained in the String Gas cosmology in which the 
emergent universe was achieved in the Hagedorn phase 
of a thermal system composed of a gas of superstrings 
E3 [35J [35] . It can be implemented by a model dubbed 
as Galilean Genesis as well |37j . Phenomenological study 
on the causal generation of primordial field fluctuations 
in this scenario was performed via the so-call conformal 
cosmology [38-43J, and the pseudo-conformal cosmology 
[44JI45]. Later, the issue of successfully transferring scale- 
invariant primordial field fluctuations to curvature per- 
turbations was discussed in [45] and [47], respectively. 
Moreover, there is a modified version of emergent uni- 
verse in which the universe has experienced a process of 
slow contraction [37] 145^50 or slow expansion [5TM56| . 
Recently, it has been also shown that emergent-type uni- 
verses, which start out with a very small Hubble rate, 
are quite generally preferred over inflationary models in 
a landscape, as described in [571 [55] , 
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In this paper, we are going to show that within the 
frame of the 4-dimensional Friedmann- Robertson- Walker 
(FEW) universe governed by standard Einstein gravity, 
the Quintom-like matter field is needed to realize a realis- 
tic model of emergent universe, which drives the universe 
from a quasi-Minkowski phase to a normal thermal ex- 
panding history. We start with a brief examination on 
this necessary condition. 

Consider a universe initially emergences from a 
nonzero minimal size and experiences enough long period 
of quasi-Minkowski phase. In this phase, the scale fac- 
tor a(t) is almost a constant and its time derivative a(t) 
is nearly zero. In order to make the universe exit this 
phase gradually, we need its second order time derivative 
ii(t) to be a very small positive value. Thus, if we use the 
Hubble parameter H (= a /a) to characterize dynamics of 
this phase, we find H — > + and H > 0. Having assumed 
Einstein gravity is still available to describe the gravity 
sector during this period, one can immediately read there 
is an effective EoS of the universe which satisfies, 



2H 

W = - 1 ~3IP << - 1 



(1) 



After exiting the quasi-Minkowski evolution gracefully, 
the universe needs to enter into the normal thermal ex- 
panding phase which requires the EoS of the whole uni- 
verse to be roughly equal to 1/3, and —1 along with the 
observable history, respectively. Therefore, this requires 
a transit of the background EoS from w < — 1 to u; > —1, 
which is exactly the Quintom behavior. 

In this paper we study particular realizations of the 
emergent universe picture by using several explicit Quin- 
tom models. The paper is organized as follows. In Sec- 
tion II, we study the general requirements of the emer- 
gent universe picture. Then in Section III, we present the 
analytic and numerical solution of the emergent universe 
from a toy model of a phenomenological Quintom-like 
fluid with a parameterized EoS across the cosmological 
constant boundary. For a much concrete model building, 
we consider a nonconventional spinor field |60j in Sec- 
tion III. Particularly we make use of an ansatz of the 
EoS and reconstruct the potential of the spinor Quintom 
analytically and numerically. Section V is devoted to a 
summary of the paper. 



II. THE REQUIREMENTS OF THE 
EMERGENT UNIVERSE 



To start, we consider a spatially flat FRW universe, of 
which the metric is given by 

ds 2 = dt 2 - a 2 (t)dx ■ dx . (2) 

Without modifying General Relativity, it is well known 
that the background dynamics in this frame follow the 
following two equations of motion, 

p + P 



H 1 



P 

3M p 2 ' 



H 



2Mj ' 



(3) 



in which M p = 1/\Z8ttG is the reduced Planck mass. In 
addition, p and P are the energy density and the pressure 
of the matter fields filled in the universe, respectively 

The emergent universe is a scenario of non-singular 
cosmology. Differing from the bouncing or cyclic modes 
where the expanding universe was preceded by a con- 
tracting phase, the picture of emergent universe requires 
the universe expands forever beginning with a finite scale 
factor in the infinite past. This implies the Hubble rate 
H cannot be negative. 

To be non-singular requires the spacetime of the uni- 
verse is geodesically complete, i.e., the affine parameters 
of geodesies are divergent in the limit of infinite past [55] • 
The null geodesic obeys the equations 



dk 



and 



dx ' = ° 



dfivk k — , 



(4) 



(5) 



where A is the affine parameter and = dx^/dX is the 
vector tangential to the geodesic. In the spatially flat 
FRW universe with the line element one can show 
that 



dk° 

~dT 



Hk° = 



(6) 



By making use of H = dlna/dt, the above equation 
yields 



k° = — oc — , dX cx a(t)dt 
dX a 



(7) 



and hence, the requirement of the completion of null 
geodesies can be expressed as 



t'=t 



a(t')dt' 



a(t'=i) 



a(t' — — oo) 



H^da] = oo . (8) 



Because the scale factor a is finite and non-negative, the 
divergence of the second integral needs Pl~ y should be 
singular at certain moment. For example, in bouncing 
cosmology, H^ 1 = oo at the bouncing point which is 
located at some intermediate time. However, for the 
emergent universe scenario, it is natural to conclude that 
H^ 1 oo when t' — > — oo, which means, the spacetime 
of the universe approaches Minkowskian in the infinite 
past. 

For time-like geodesies, we can think that they are 
worldlines of some free particles with masses. Consider 
a wordlinc of a free particle with mass m, the proper 
time s itself is the affine parameter of the wordlinc. The 
four-momentum is defined as 



m 



dx^ 
ds 



(9) 
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which obeys g^P^P" = to 2 . The condition for the 
geodesic completion is obtained from the geodesic equa- 
tion, which is 




ff-v + ^r 1/2 <foi = > (10) 



where C p — PiPi is a non- negative constant. This 
also requires the universe starts at infinite past from a 
Minkowski spacetime where H^ 1 = 00 or H = 0. The 
evolution of a Minkowski spacetime to an expanding uni- 
verse requires the H increases for some time during which 
H > 0, this implies the equation of state w < — 1 and 
the null energy condition was violated. At later time the 
universe should enter into the radiation dominated phase 
in which w — 1/3. So w crosses —1 at some intermediate 
time. That is to say, the matter dominating the universe 
has Quintom behavior. 

Furthermore, in the emergent universe the Hubble rate 
reached a maximum H max at the crossing point. It de- 
fines a mass scale the universe cannot go beyond. If H max 
is much smaller than the Planck mass, the validity of 
the general relativity as a low energy effective theory is 
guaranteed and the corrections of quantum gravity are 
suppressed. 
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FIG. 1: Plot of the evolutions of the scale factor a, the Hubble 
parameter H and the background EoS w as a function of 
cosmic time in the case of Quintom fluid. In the numerical 
calculation, we take a = 3 and C = 10. All dimensional 
parameters are of Planck units. 



III. EMERGENT UNIVERSE WITH A 
QUINTOM-LIKE FLUID 

As a first example, we illustrate the possibility of ob- 
taining the emergent universe solution in a phenomeno- 
logical Quintom fluid described by the following EoS: 



w(t) = i - ^e-^* 
w 3 3 



(11) 



where a is a positive-valued dimensionless parameter. 
In this particular parametrization, we have assumed the 
universe can automatically enters the radiation domi- 
nated period after the quasi Minkowski expansion. Thus, 
one can see the EoS approaches 1 /3 when the cosmic time 
t goes to positive infinity. Moreover, when t approaches 
to a far past moment, the EoS would have fallen down to 
negative infinity very soon due to the expression of the 
exponential term. 



To substitute the parametrization of the EoS (111 into 



the Friedmann equations, one can solve out the explicit 
solution to the Hubble parameter as follows, 



H(t) 



M p e 



aM„t 



1 + (C + 2M p t)e 



(12) 



where C is a constant which can determine the energy 
scale of the Hubble parameter at the occurrence of emer- 
gent universe. Further, one obtains the energy density 
and the pressure of the Quintom fluid governing the uni- 
verse as well. As a consequence, the evolution of the scale 
factor in this model can be numerically integrated out 



as shown in Fig. [T] From the figure, one can explicitly 
read that the scale factor a approaches to a non- vanishing 
minimal value in the limit of far past, and connects the 
thermal expansion when t — > +00. 

Notice that, the occurrence moment ts of the emergent 
universe scenario can be characterized by the moment of 
w crossing —1. Thus, by solving w(tE) = — 1, one gets 



t E 



1 



aM, 



log(") 



(13) 



which in our explicit example takes the value of 0.14 ap- 
proximately. At this moment, the Hubble parameter ar- 
rives at the maximal value 



H E = 



M„ 



C+f(l-log(f)) 
IM P . The 

Einstein gravity in this case. 



(14) 



which is around 0.1M p . Therefore, we are able to trust 



IV. EMERGENT UNIVERSE AND SPINOR 
QUINTOM 

In this section, we consider a class of Quintom model 
described by a nonconventional spinor field. To begin 
with, we simply review the background dynamics of a 
spinor field which is minimally coupled to Einstein's grav- 
ity (see Refs. [611 EZ] for detailed introduction and see 
[BniISS, 64 for recent phenomenological study in cosmol- 
ogy)- 
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A. Algebra of a cosmological spinor 

Following the general covariance principle, a connec- 
tion between the metric g^ and the vierbein is given 
by 



9^v^a e b — Vab 



(15) 



where denotes the vierbein, g^ is the space-time 
metric, and r/ab is the Minkowski metric with r/at — 
diag(+l, — 1, — 1, — 1). Note that the Latin indices rep- 
resents the local inertial frame and the Greek indices 
represents the space-time frame. 

We choose the Dirac-Pauli representation as 



7° = 



1 
-1 



7 = 



0* 

-at 



(16) 



where cr, is Pauli matrices. One can see that the 4x47° 
satisfy the Clifford algebra {7°,7 fc } = 2r) ab . The 7 Q and 
provide the definition of a new set of Gamma matrices 



<1° 



(17) 



which satisfy the algebra {r M ,T' / } = 2g lxv . The genera- 
tors of the Spinor representation of the Lorentz group can 
be written as S ab = \[y a , 7 b ]. So the covariant derivative 
are given by 



Dftip = d^ip - ipCl^ , 
where the Dirac adjoint is defined as ip 



0+7°. 



(18) 
(19) 

The 

4x4 matrix f2 M = ^u!f ia b^ ab is the spin connection, where 
Uf+ab = e V aS t± e vb are Ricci spin coefficients. 

By the aid of the above algebra we can write down the 
following Dirac action in a curved space-time background 



d 4 x e {-(^D^ - D^T^) - V$iP)] .(20) 



e is the determinant of the vierbein e" and V stands 



Here 

for the potential of the spinor field ip and its adjoint ip. 
Due to the requirement of covariance, the potential V 
only depends on the scalar bilinear ipip and "pseudo- 
scalar" term ip^ip. For simplicity we drop the latter 
term and only assume V = V(ipip). 

Varying the action with respect to the vierbein we 
obtain the energy-momentum-tensor, 



T, 



C/xa 5 Sip 

-g^C-ip ■ (21) 

On the other hand, varying the action with respect to 
the field ip, if> respectively yields the following equations 
of motion, 



iT^D^ip - V'ip 
iD^T" + V'ip 



. 
. 



(22) 
(23) 



where V = dV/d(ipip) denotes the derivative of the 
spinor potential with respect to iptp. 

We deal with the homogeneous and isotropic FRW 
metric. Correspondingly, the vierbein are given by 



1 



(24) 



Assuming the spinor field is space-independent, the equa- 
tion of motion reads i^°(ip+^Hip) — V'ip = 0, where a dot 
denotes a derivative with respect to the cosmic time and 
H is the Hubble parameter. Taking a further derivative, 
we can obtain: 



ipip 



N 



(25) 



where N is a positive time-independent constant. 

From the expression of the energy-momentum tensor 
in Eq. (21 ), one can obtain the expressions of the energy 



density and the pressure of the spinor field as follows, 



Pn, = r u = v , 

Pi, = -T? = V'i^iP - V 



(26) 
(27) 



As a consequence, the EoS of the spinor field is given by 

V'ipip 



Pi, 

Wil, = — 

Pi, 



-1 



V 



(28) 



B. Reconstruction of Spinor Quintom realizing 
emergent universe scenario 

The above formulae show that a cosmological spinor 
field might realize its EoS to cross the cosmological con- 
stant boundary when the sign of V changes. For a con- 
ventional spinor with its potential taking the form of 
mtpip, its EoS is exactly zero which coincides with that 
of normal non-relativistic dust matter. Thus, one has to 
consider a nonconventional form of the potential for the 
cosmological spinor to achieve the Quintom scenario. 

Interestingly, observing the form of the Friedmann 
equations and the formulae of the spinor field, we can 
see that a model of spinor Quintom can be reconstructed 
to fulfill the scenario of emergent universe. Suppose a 
fixed form of the EoS such as what we have illustrated in 
the case of Quintom fluid. We then read the expression 
of the Hubble parameter H and its time derivative H, 
and thus can further derive the evolutions of V and V 
as functions of cosmic time. By making use of the rela- 
tion that ipip oc a -3 , the form of V can be derived out 
explicitly. We will do the reconstruction in the following 
context. 

Note that, a realistic cosmological evolution requires 
a thermal expansion after the primordial period. Af- 
ter radiation dominated phase, the universe enters into 
a matter dominated era. However, we usually do not ex- 
pect that a fermion field could be responsible for a large 
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amount of radiation 1 . Therefore, for a simple and natu- 
ral choice, we parameterize the form of EoS for the spinor 
field as follows, 



2a 



-aM„t 



(29) 



To observe this form, one can see that the EoS falls 
into negative infinity when t <C — 1, but approaches 
at t — > +oo. This parametrization can give rise to a 
emergent universe solution with a dust-like expansion fol- 
lowing the quasi-Minkowski phase without the radiation 
domination. 



Substituting Eq. (29 1 into the Friedmann equations, 



one can solve out the Hubble parameter as a function of 
cosmic time: 



H(t) 



M n e a 



1 + (C+ fM P t)e« M p* 



and thus the potential of the spinor evolves as 



V 



12M*e 



4„2aAf„t 



(2 + 2(C + 3Mpt)e QM p*) 2 



(30) 



(31) 



Similar to the case of Quintom fluid, the coefficient C is 
a integral constant which is used to determine the energy 
scale of the occurrence of emergent universe scenario. 

By numerically solving the Friedmann equations, one 
can solve out the evolution of the scale factor along cos- 
mic time. We show the evolutions of the scale factor a, 
the Hubble parameter H and the EoS w in Fig. [2j 

After having known the evolution of the scale factor, 
one can numerically obtain that of the scalar bilinear ipip 
due to the relation that 'iptp oc a~ 3 . In addition, the 
energy density of the spinor field only depends on the 
potential V and thus one can derive V evolving along 
cosmic time. They are shown in Fig. |3[ From the upper 
panel of this figure, one can see that ipip is a monoton- 
ically decreasing function and approaches a constant in 
the emergent universe period. Note that, we use log scale 
to show the wide range of scales for V in the longitudinal 
coordinate of the lower panel of Fig. |3j where it can be 
seen that V evolves as an exponential form in the phase 
of quasi-Minkowski. 

Since the evolutions of tptj) and V were already ob- 
tained in above numerics, we can further derive V as a 
function of x/ji/j. This is achievable as ipi/i is decreasing 
monotonically along the cosmic time and then can lead 
to a inverse function t = t(ip?p). To combine t = t(?pip) 
and V = V(t), we then numerically solve out V = V("ijjijj) 
as shown in Fig. [4] 

Interestingly, from Fig. [4] one can find that V is a 
linear function of tpip at the regime of small values. 



1 During radiation dominated phase, the main contribution comes 
from the gauge photons, and only a few part is from nearly mass- 
less hot neutrinos. 




FIG. 2: Plot of the evolutions of the scale factor a, the Hub- 
ble parameter H and the background EoS w as functions of 
cosmic time in the case of spinor Quintom. In the numerical 
calculation, we take a = 3 and C = 10 which are the same as 
as in the case of Quintom fluid (Fig. [IJ, and N = 100 for the 
scalar bilinear ipip. All dimensional parameters are of Planck 
units. 



This implies that the spinor Quintom recovers the con- 
ventional form of a massive fermion with its potential 
V ~ mipip at low energy limit. Moreover, in the UV 
limit, one can derive the asymptotical form of the scale 
factor as follows, 



a(t)^a E (l + Ce aM ^)^c 



(32) 



and thus the asymptotical forms of iptp and V. To take 
the inverse function of ipip(t), eventually we can get the 
approximate expression of the potential in the phase of 
quasi-Minkowski: 



3M p 4 
-(jT 



2E ( lv } 



(33) 



where being the minimal value of the scale factor in 
the emergent universe scenario. 

Note that, the potential for the spinor Quintom in 
low energy limit is very normal in quantum field theory, 
which is merely a mass term. However, the potential 
becomes very nontrivial at the high energy scale, which 
takes the form (33). Although this form is purely phe- 



nomenologically constructed, one can find that it implies 
a condensate of the tachyonic spinor field in UV regime 
and thus is expected to be related to the spinor formalism 
in open string field theory [551 166j . 
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FIG. 3: Plot of the evolutions of the scalar bilinear ijjip and 
the potential V as functions of cosmic time in the case of 
spinor Quintom. In the numerical calculation, we take the 
values of a, C and N the same as in Fig. [2] All dimensional 
parameters are of Planck units. 
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FIG. 4: Plot of the potential V as a function of •tpip in the case 
of spinor Quintom. In the numerical calculation, we take the 
values of a, C and N the same as in Fig. [2] All dimensional 
parameters are of Planck units. 



V. CONCLUSION 

As well known, the conception of Quintom scenario was 
originally from the study of dark energy physics, which 
shows the EoS of dark energy might cross the cosmo- 
logical constant boundary |14j . This scenario is mildly 
favored by a group of cosmological observations [57J |5S] , 
and its dynamics were extensively analyzed in a number 
of works (for example see |69l - f76"] for phenomenological 
study and see [19 for a comprehensive review). How- 
ever, it was soon found that the Quintom behavior has 
many significant implications to early universe physics. 
Namely, it can give rise to nonsingular bouncing and 
cyclic [77H75] solutions when applied to high energy scale. 

To conclude, we in the current paper have studied the 



realization of emergent universe scenarios in the presence 
of Quintom matters. In the literature there have been a 
lot of efforts in building the models of emergent universe 
to avoid the big bang singularity. Also there are efforts 
on investigating the perturbation analysis in the phase of 
emergent universe assuming the gravity sector is still de- 
scribed by standard Einstein gravity. However, we in the 
present paper point out that a model giving rise to the 
emergent universe scenario has to be of Quintom behav- 
ior within the standard 4-dimensional FRW framework 
under the assumption of Einstein gravity. 

In explicit realizations, we first considered the Quin- 
tom fluid with a parameterized form of EoS to illustrate 
its possibility. After that, we have studied in detail the 
model of spinor Quintom. As the model of spinor Quin- 
tom possesses a very nice algebra relation in the curved 
spacetime, we are able to reconstruct a suitably chosen 
form of the potential and force it to give rise to the emer- 
gent universe scenario. Moreover, this type of model 
is able to exit the phase of quasi-Minkowski expansion 
gracefully and enter a normal expanding phase domi- 
nated by dust matter. 

The scenario proposed in the current letter and its im- 
plementation in the detailed model has some remarkable 
properties in phenomenological applications. 

First of all, it is interesting to realize the Quintom sce- 
nario in other effective field models, such as the Galileon- 
type fields. In the original model of Galilean Genesis, it 
was found that the universe could hardly exit the emer- 
gent state smoothly since its equation of state was unable 
to cross —1 from below to above. According to our anal- 
ysis, one can expect a much improved model of Galilean 
Genesis which can take the advantage of Quintom sce- 
nario to exit the primordial era smoothly. 

Additionally, a crucial issue in the cosmology of emer- 
gent universe is the processing of cosmic perturbations 
throughout the quasi-Minkowski expanding phase. In the 
literature, there have been extensive studies on the gen- 
eration of primordial power spectrum. As a first step, 
we in the present paper only considered the background 
evolution but requires the energy scale of the emergent 
universe to be much lower than the Planck scale, so that 
the assumption of Einstein gravity is trustable in our 
model. However, we should be aware of that since a 
cosmic spinor was introduced to realize the scenario of 
emergent universe, it may seed some unwanted fluctua- 
tions modes such as vector modes through gauge inter- 
actions. Such a complete perturbation analysis of the 
emergent universe Quintom model lies beyond the scope 
of the present work and it is left for future investigation. 
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pressure of the spinor field: 

= ~(W>Do1>-Dofrfil>)-V(W>) 

= V'^ip - V , (37) 

where we have used the assumption diip = in the second 



line and applied Eqs. (34 1 and (35) in the last line. 



B. The stability issue 



Appendix 

In the first part of this Appendix, we derive the explicit 
expressions for the energy density and pressure provided 
in Eqs. (26) and (27), respectively. In the second part of 



the Appendix, we address on the stability of our model 
throughout the cosmological evolution. 



A. Energy density and pressure 

We assume the spinor field is space-independent but 
evolves along cosmic time. Then the equations of motion 



Regarding the stability issue of a cosmological model, 
there are mainly two possible concerns. One is to check 
whether or not the model leads to a ghost degree of free- 
dom. Obviously, this issue does not exist in our emergent 
universe model which is realized by a spinor Quintom. 
As one can see from the action, the kinetic term for the 
spinor field is very standard without any modification. 
Thus, there does not exist any extra degree of freedom in 
our model. The other is to study the evolution of cosmo- 
logical perturbations seeded by the spinor field and make 
sure the backreaction of perturbations to the background 
is controllable. 

Here we would like to show the perturbation theory of 
Spinor Quintom crudely. In order to simplify the deriva- 
tive, we would like to redefine the spinor as ipN = a^ip, 
and perturb it as 



given in (22) and (23) are simplified as follows, 



(38) 



il°(^ + ^H^) -V'tp = 
i($ + |-ffV07° + V'tjj = 



(34) 
(35) 



According to the expression of the energy-momentum 
tensor given in Eq. (21), we calculate the expression of 



the energy density of the spinor field: 



-L n 



2a 



(36) 



where in the last equality we used the assumption of scale 
independence di%p — 0. Furthermore, we can compute the 



around the homogeneous background. Then perturbing 
the equation of motion of the spinor field, one can obtain 
the perturbation equation as follows, 



d 1 

dr 2 



6ip N - V 2 5ip N + a 2 [V' 2 
= -2a 2 V'V"S(^i/j)i/} N 



vf(HV - 3HV"i}tp)]Sxj; N 
-i^d^aV"5{^)]^ N , (39) 



where r is the conformal time defined by dr = dt/a. 
Since the right hand side of the equation decays propor- 
tional to a~ 3 or even faster, we can neglect those terms 
throughout the evolution of the universe for simplicity. 

From the perturbation equation above, we can read 
that the sound speed is equal to 1 which eliminates the 
instability of the system in short wavelength. Moreover, 
when the equation of state w crosses —1, we have V — 
at that moment and the eigen function of the solution to 
Eq. ( 39 ) in momentum space is a Hankcl function with an 



index ^. Therefore, the perturbations of the spinor field 
oscillate inside the hubble radius. This is an interesting 
result, because in this way we might be able to establish 
the quantum theory of the spinor perturbations, just as 
what is done in inflation theory. 
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Moreover, in the emergent universe phase one can in- 
sert the expressions of scale factor ( 32 ) and Hubble pa- 



rameter (30), and apply the approximate form of the po- 
tential ( 33 ) into the perturbation equation ( 39 1 . Then 
one can get H ~ e QMpt and V - • " u ' 
ponentially suppressed but V" ■ 



e "*«f>« which are ex- 
constant. This implies 



that the perturbations of the spinor field are almost vac- 
uum fluctuations in the phase of emergent universe which 
take the form of plane wave function approximately. As a 
consequence, the backreaction of perturbations are neg- 
ligible in our model. 
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